Application of finite element method to two-dimensional multi-group neutron transport equation in cylindrical geometry

藤村 統一郎; 筒井 恒夫; 堀上 邦彦; 中原 康明; 大西 忠博*

Journal of Nuclear Science and Technology, 14(8), p.541 - 550, 1977/08

https://doi.org/10.3327/jnst.14.541

先に、二次元(r,z)円柱体系における多群中性子輸送問題を有限要素法で解くアルゴリズムが開発され、簡単なモデルによる計算もなされた(日本原子力学会、昭和49年秋の分科会、同昭和50年年会での口頭発表)。 有限要素法は(r、z)面内の正規長方形小領域上の空間変数に応用されている。 本稿では、そのアルゴリズムのうち、双一次または双二次多項式を基底として用いた不連続法およびその計算結果について述べる。 原子炉の現実的な体系を中心としたいくつかの数値例が示されるが、双二次近似による解は精度も良く、粗いメッシュのときでも数値的に安定である。 また、汎用的なダイヤモンド差分法によるコードとの比較もなされ、また数値計算の結果を通じて不連続法の利点が示されている。